Multimode Memories in Atomic Ensembles 
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The ability to store multiple optical modes in a quantum memory allows for increased efficiency 
of quantum communication and computation. Here we compute the multimode capacity of a variety 
of quantum memory protocols based on light storage in ensembles of atoms. We find that adding a 
controlled inhomogeneous broadening improves this capacity significantly. 
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Quantum memories with the ability to store multiple 
optical modes [l|, [3] offer advantages in speed and robust- 
ness over single mode memories when incorporated into 
quantum communication [3, |j, Q and computation [a Q 
schemes. For example, dual rail encoding 8] and photon 
loss codes [9| can be integrated into ensemble based com- 
putation schemes if multimode storage is available. It is 
therefore useful to characterize the multimode capacity 
of a given quantum memory protocol; this is the problem 
we address in this paper. 

Any quantum memory must operate as a linear map, so 
that quantum superpositions are maintained. Therefore 
analysis and optimization of the storage efficiency and 
multimode capacity of any quantum memory can be un- 
derstood using the well-developed tools of linear algebra. 
Of central importance is the singular value decomposi- 
tion (SVD), or Schmidt decomposition [lO|, [ill, which 



exists for all linear maps. In this paper we introduce a 
universal method for calculating the multimode capacity 
of a quantum memory, based on the SVD of the Green's 
function for the memory interaction. We apply this tech- 
nique to memories with no inhomogeneous broadening, 
and we find that their multimode capacity scales poorly 
with the optical depth of the ensemble. We go on to 
calculate the multimode capacity of protocols based on 
controlled reversible inhomogeneous broadening (CRIB) 
[121 Il3| , and find that their capacities scale more favor- 
ably. We then consider a modification to the well-known 
Raman memory protocol IJ, llSl . Il6l llTI ] , in which a lon- 
gitudinal inhomogeneous broadening is applied 18|, ll9| ; 
we show that this modified protocol exhibits the same 
favorable multimode scaling as CRIB. Finally we exam- 
ine the scaling of the recently proposed atomic frequency 
comb (AFC) protocol [l|], which was designed with mul- 
timode storage in mind; its performance is the best out 
of the protocols we analyzed. 

In the following we deal exclusively with quantum 
memories based on ensembles of atoms, or atom-like ab- 
sorbers, although we neglect memories based on mea- 
surement and feedback 20|, |2l|, [23, l23|, \2A\- In a typical 
ensemble quantum memory, an incident signal field with 
temporal envelope Ai,^{t) is converted to a stationary 
excitation known as a spin wave, which is distributed as 



a function of position z within the ensemble with ampli- 
tude B{z) (see Figure [T] (a)). The quantum memory map 
takes the form 
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where the kernel K is a, Green's function for the storage 
interaction. In some cases, the form of K can be derived 
in closed form. The general strategy is to solve the lin- 
earized Maxwell-Bloch equations describing the coupling 
of the signal to the atoms — the Fourier transform is an 
invaluable tool in this connection. Even when an analytic 
solution is not possible, K can always be constructed nu- 
merically. In fact, although we use a one dimensional 
treatment here and in what follows, the same arguments 
apply if transverse coordinates are included, although the 
numerical problem becomes significantly larger in this 
case. Discussions of multimode storage using the trans- 
verse spatial degrees of freedom can be found elsewhere 
13) la, l25|. Using the SVD, we can always express the 
Green's function as a convex sum 
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where the sets of modes {ipk}, {'t'k} each form orthonor- 
mal bases, the former for the space of spin wave am- 
plitudes — the 'output space' — and the latter for the 
space of signal envelopes — the 'input space'. Note that 
the SVD does not require the kernel K to be Hermitian: 
in general, losses can make the kernel complex and non- 
Hermitian. The positive real singular values Xk couple 
each input mode <pk to the corresponding output mode 
i/jk- The probability that a single photon in mode 0^ is 
stored by the memory is A^ . The optimal performance of 
a memory is achieved when the storage and retrieval pro- 
cesses are simply time reverses of each other [15|, l26|, l27| , 
in which case the retrieval interaction is described by the 
conjugated Green's function K*. The retrieval efficiency 
for the fc*'' mode is then the same as the storage effi- 
ciency A^., and so the total efficiency of storage followed 
by retrieval — the 'memory efficiency' — is given by A^. 
The maximum number of modes of the optical field 
that can be stored by the memory is just the number 



of non-negligible singular values. We note that the or- 
thonormal basis {(/>fc} is the 'natural' input basis for the 
memory interaction, as determined by its dynamics. It 
represents the best possible encoding, and the multi- 
mode capacity can only be reduced by choosing some 
other encoding, such as time or frequency bin encoding, 
for the input modes. To quantify the multimode capac- 
ity of a memory, we label the Xk in order of descending 
magnitude, and define the average memory efficiency for 
the first k modes, A^ = J2j=i X'j/k. Let at be equal 
to 1 when A^ exceeds some threshold efficiency 9, with 
ttfe = otherwise. Then we define the multimode capac- 
ity N = '^f.ctk as the number of modes with average 
memory efficiency greater than 6. In this paper we use 
6 = 70% for all calculations; this choice makes the mul- 
timode scaling apparent for experimentally achievable, 
and numerically tractable, parameters. We have checked 
that the scaling of N does not depend on the choice made 
for 6, although imposing higher thresholds reduces N, as 
might be expected. 
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FIG. 1: (a): An incident signal field is mapped to a stationary 
spin wave inside an ensemble quantum memory, (b): The 
multimode scaling of unbroadened ensemble memories (solid 
line) alongside the tCRIB (dashed line) and ICRIB (dotted 
line) protocols, with a threshold 9 = 70%. 



The SVD of K is related to the spectral decomposition 
of the normally and anti- normally ordered products {c.f. 
Theorem 5.4 in ^) 



Kn{t,t') = 



K*{z,T)K{z,T')dz, 
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K{z,T)K*{z',T)dT. 



(3) 



Direct substitution of Eq. ^ into the above expressions 
shows that the eigenvalues of both Ka and K^f are given 
by fik — A^; the capacity N can therefore be directly 
calculated from these eigenvalues. Here we normalized 
the longitudinal coordinate so that z runs from to 1 . In 
some cases we may have access to an explicit form for the 
Green's function Kt for the entire memory interaction, 
that is, storage followed by read out. Suppose that we 
have 
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where Ao^t is the temporal profile of the signal field re- 
trieved from the memory. This is, in some sense, more 



complete than the storage Green's function in Eq. ([T]), 
since it tells us directly about the connection between 
the modes that can be stored and the modes that can 
be retrieved. The SVD can be applied to the bivariate 
function Kt precisely as it was to K in Eq. ^, and the 
singular values i>k extracted. In the optimal case that re- 
trieval is the time reverse of storage, we have Kn = Kt, 
and so for consistency we define the multimode capacity 
in terms of the Vk by making the replacement A^ —>■ Vk- 

Finally, we note that the singular values are invari- 
ant under unitary transformations. In particular, we are 
free to Fourier transform position to spatial frequency 
(z — > k) or time to angular frequency (r — > w); the sin- 
gular values of a kernel K are unchanged by such trans- 
formations, and so the multimode capacity can equally 
well be calculated from the transformed kernel K. 

It has previously been shown [17l] that optimal storage 
in an ensemble memory without inhomogeneous broad- 
ening — whether b y using electromagnetically induced 
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storage, or simply linear absorption 

by the kernel 



or Raman [3, [H, 0, [Tl | 
is characterized 
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where d is the resonant optical depth of the ensemble 
17 1 , and Iq is a modified Bessel function [34I . Using this 
result, we calculate the multimode capacity as a function 
of the optical depth d (see Figure [T](b)). A numerical fit 
reveals a square-root scaling N '^ ^/d/3. We can iden- 
tify this as the scaling of the absorption linewidth of an 
ensemble of 2-level atoms, each with a Lorentzian line- 
shape. Since large optical depths can be challenging to 
achieve in the lab, this scaling makes unbroadened en- 
semble memories unattractive candidates for multimode 
storage. Can this poor scaling be improved? Clearly we 
can always store N modes if we have A^ separate memo- 
ries, as long as each has sufficient optical depth dg such 
that Af > 9. Then the total optical depth is simply 
d = Ndg, so that N scales linearly with d. Below we 
show that in fact separate ensembles are not required for 
linear scaling. A single physical memory can be made 
'parallel' by redistributing atoms in frequency. 

We now calculate the multimode capacity for memory 
protocols based on controlled reversible inhomogeneous 
broadening (CRIB) [Ij, [l3|, in which a spatially vary- 
ing external field is applied to an ensemble in order to 
artificially broaden the range of frequencies absorbed by 
the atoms. A signal field is absorbed resonantly by the 
broadened ensemble. Some time later, the polarity of 
the external field is reversed, so that the inhomogeneous 
broadening profile is flipped, with red-detuned atoms be- 
coming blue-detuned and vice versa. The optical dipoles 
eventually re-phase, and the signal field is re-emitted. 
The optimal efficiency is achieved if the signal field is re- 
emitted in the backward direction, but phasematching 



considerations sometimes favor forward emission. In the 
following we consider ideal backward emission, since for- 
ward emission cannot perform any better. The artificial 
broadening can be applied in a direction transverse to 
the propagation direction of the signal field to be stored 
(tCRIB), or longitudinally — parallel to the signal beam 
(ICRIB). As in the rest of this paper, we use a one dimen- 
sional propagation model. Since tCRIB is laterally asym- 
metric, coupling between the transverse spatial structure 
and the spectrum of the signal is potentially important. 
Nonetheless the results presented below represent an up- 
per bound on the performance of tCRIB. The Green's 
function relating the input and output spectra of the sig- 
nal field is given by [13, l30| , 
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where the 
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lineshape function is /(^) = 

..^ s dA, with p(A)dA the proportion of 
atoms whose optical resonance is shifted by a detuning 
A away from its unbroadened frequency. Here and in 
what follows 27 is the homogeneous linewidth of the 
unbroadened transition. Note that d is the total optical 
depth of the ensemble — if the broadening were removed, 
the total resonant optical depth of the ensemble would 
be given by d; this definition enables a direct comparison 
with unbroadened memory protocols. The capacity of 
tCRIB is plotted in Fig. ([1]) (b), optimized with respect 
to the spectral width Aq of the applied broadening, 
assuming a rectangular broadening profile. The scaling 
is manifestly linear, with N ~ d/2h. We found that the 
optimal width also scales linearly, Aq^ /^ ~ 9(i/5. These 
results can be explained qualitatively by considering 
spectral modes of bandwidth 7. Every such mode 
requires a resonant optical depth dg > 10. To store 
N modes 'side by side' in frequency we should have a 
broadening Ao ~ N^ and a total optical depth d ~ Ndg, 
so that N rises linearly with d, provided we increase Aq 
at the same time. Next we consider ICRIB. 

The storage kernel for ICRIB, where the atomic res- 
onance frequency varies linearly with the longitudinal 
coordinate z, can be calculated in terms of the spatial 
frequency k of the spin wave. In this case the spin wave 
represents the atomic polarization. We obtain [33| 
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where Aq is the total spectral width of the broadening 
applied, and (3 = dj/Ag. Note that k is dimensionless, 
because of the normalization of z. Unfortunately this 
kernel is singular when k — —Aqt, so a numerical SVD 
fails. We can, however, directly construct the Green's 
function in terms of z and t, which is not divergent, by 
solving the optical Maxwell-Bloch equations numerically. 
In Fig. ID) (b) we plot the multimode capacity, optimized 



over the broadening width Aq , found from such a numeri- 
cally constructed Green's function. We used spectral col- 
location [3l| for the spatial propagation, and a second- 
order Rungc-Kutta method for the time stepping. The 
capacities of tCRIB and ICRIB are nearly identical, and 
our numerics suggest that the multimode scaling of both 
these protocols is the same. 

Having shown that artificially broadened memories ex- 
hibit superior multimode scaling, we now consider a mod- 
ification to the standard Raman protocol, in which a lon- 
gitudinal broadening is applied to the storage state. The 
hope is that such a protocol will combine the advantages 
of Raman storage, namely broadband, tunable operation, 
insensitive to unwanted inhomogeneous broadening, with 
the superior multimode capacity of ICRIB. This type of 
'hybrid' storage has recently been experimentally demon- 
strated using magnetic gradients in a warm Rubidium 
vapor [l8| . The atoms comprising such a memory have a 
A- type structure 15|, [iTI, ll8| , with an excited state cou- 
pled strongly to both a ground and a storage state. A 
direct transition from the ground to the storage state 
is forbidden, and therefore storage of an incident signal 
field is mediated by an off-resonant Raman interaction 
with the excited state, the strength of which is controlled 
by the application of an intense control field coupling the 
storage and excited states. An external field is applied so 
that the energy of the storage state varies linearly along 
the length of the ensemble, covering a spectral width Aq. 
The spin wave B in this protocol represents the coher- 
ence betwe en grou nd and storage states. Under adiabatic 
conditions 15l.ll7|. the following expression approximates 



the memory map. 
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where 5(fc,T) = [d-y + iVx (fc-|-Aor)]"^ Here F = 7 + iA 
is the complex detuning, with A the common detuning of 
the signal and control fields from the excited state. Q, is 
the slowly varying Rabi frequency of the control field. 
Figure ^ (a) shows the multimode capacity for this 
broadened Raman protocol, operated with a Gaussian 
control field. Changing the temporal profile of the con- 
trol pulse, while maintaining a constant energy, affects 
the shapes of the input modes {0^}, but not the multi- 
mode capacity. As anticipated, the square-root scaling 
of the unbroadened memory — smaller than predicted 
by Eq. ^ because we operate far from resonance — is 
transformed into the linear scaling of the CRIB protocols 
by apphcation of a broadening. We found N ^ d/300 
with Aq'' /7 ~ d/77; in general larger optical depths are 
required to achieve the same capacity for the Raman pro- 
tocol, again due to the large detuning from resonance. 

Finally we turn to the recently proposed atomic fre- 
quency comb (AFC) protocol [1[. This protocol is sim- 
ilar to tCRIB, except that the inhomogeneous profile 
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FIG. 2: (a): The multimode scaling for the Raman protocol 
with an applied broadening (solid line) and without (dashed 
line). Here we set f7(r) = Vwd-je-^^^''^''^ and A = ^9057. 
(b): The capacity of AFC with several different numbers of 
comb teeth, indicated by the numbers in the plot. We used a 
threshold 6 = 70%. 



p(A) takes the form of a comb with spectral width 
Aq, comprised of M equally spaced resonances, each 
with optical depth d: p(A) — J2i=iH'^ ~ ^j)^ with 
Sj = -Ao/2+ (j - l)Ao/(M- 1). The frequency comb is 
prepared from an initially uniform naturally inhomoge- 
neously broadened ensemble by optical pumping, which 
essentially removes atoms that would absorb at frequen- 
cies between the comb teeth. More teeth can be added 
by pumping out fewer atoms, which does not require any 
change in the size or density of the ensemble. Therefore 
the important physical resource is the depth d associated 
with each tooth, since this is set by the ensemble density. 
Note, however, that the total optical depth dtot = Md, 
found by summing the contributions from all the teeth, 
increases as we add more teeth. This explains the re- 
markable multimode scaling discussed below. The inho- 
mogeneous broadening does not need to be 'flipped' be- 
fore retrieval, as it does for CRIB, since its discrete struc- 
ture leads to periodic rephasing of the atomic dipoles. 
The Green's function for storage, followed by backward 
retrieval, using the same notation as in Eq. (O, is [ll.l3g| 
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In Fig. ([2]) (b) we plot the multimode capacity of AFC 
for several values of M. For every M the capacity ex- 
hibits the square-root scaling characteristic of unbroad- 
ened memories, but A'^ ex M, and can be made arbitrar- 
ily large simply by adding more teeth to the comb. This 
increases dtot, but requires no increase in the ensemble 
density. Provided d > dg '^ 10 the multimode capacity 
is, in principle, infinite. Of course, as with all of these 
protocols, N can never exceed the total number of atoms 
in the ensemble. More importantly, however, as more 
teeth are added the comb width Ao must be increased to 
ensure the teeth are well-separated; we found a finesse 
^ F = Ao/2j{M - 1) > 30 is required. N is therefore 
limited because Aq is bounded by the spectral width of 
the initial inhomogeneous line. 

In conclusion, we have developed a universal approach 



to quantifying the ability of ensemble quantum memories 
to store multiple spectral/temporal modes of an optical 
field. We showed that the multimode capacity of stan- 
dard FIT and Raman memory protocols scales with va, 
while the capacity of CRIB protocols scales linearly with 
d. We considered the capacity of a modified Raman pro- 
tocol in which a longitudinal broadening is applied to 
the storage state, and found that this linear scaling is 
reproduced. Finally we considered the AFC protocol; its 
capacity is not limited by d, provided a certain threshold 
is exceeded. 
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